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Abstract 

Expressions in terms of common reactor operating parameters were derived for the ratio of 
the Grashof number to the Reynolds number, Gr/Re, the ratio of the Grashof to the square of 
the Reynolds number, Gr/Re 2 , and the Rayleigh number, Ra. Values for these numbers were 
computed for an example horizontal CVD reactor and compared to numerical simulations to 
gauge their effectiveness as predictors of the presence or absence of transverse and longitudinal 
rolls in the reactor. Comparisons were made for both argon and hydrogen carrier gases over 
the pressure range 2-101 kPa. Reasonable agreement was achieved in most cases when using 
Gr/Re to predict the presence of transverse rolls and Ra to predict the presence of longitudinal 
rolls. The ratio Gr/Re 2 did not yield useful predictions regarding the presence of transverse 
rolls. This comparison showed that the ratio of the Grashof number to the Reynolds number, 
as well as the Rayleigh number, can be used to predict the presence or absence of transverse 
and longitudinal rolls in a horizontal CVD reactor for a given set of reactor conditions. These 
predictions are approximate, and care must be exercised when making predictions near 
transition regions. 
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1. INTRODUCTION 


Chemical vapor deposition (CVD) is an important method for depositing solid films from 
reactive gaseous precursors. It finds applications in many industries, particularly when high 
purity films are required. The flow patterns inside the reactor can have a significant effect on 
the properties of the deposited material, making an understanding of such flows crucial to the 
control of the process. Mathematical models are useful for examining how changes m the 
reactor’s operating parameters affect the flow inside the reactor. The results obtained can be 
used as a guide to experimental work aimed at attaining material with required characteristics. 
General purpose commercially available codes greatly ease the work of constructing such 

models. 

Two commonly observed types of flow patterns in horizontal CVD reactors are transverse 

rolls, which have axes perpendicular to the main flow, and longitudinal rolls, which have their 

axes parallel to the direction of the main flow [1]. Longitudinal rolls result in a spiral motion 

of the fluid as it travels down through the reactor. Transverse and longitudinal rolls within the 

reactor can compromise the uniformity of the deposited material by contributing to local 

variations in the temperature and species concentrations. In the case where such rolls exist and 

the Schmidt number, 5c = (where v is the kinematic viscosity and D i m is the 

m 

diffusivity of species i in the mixture), is much greater than one, any particles that form in the 
gas phase are more likely to grow if caught in these rolls, since their residence time in the 
reactor is increased. Eliminating such flow patterns should increase the uniformity of the flow, 
thereby contributing to increased deposition uniformity. Changes in the flow pattern can be 
accomplished by altering the reactor operating conditions. 
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Two-dimensional flow calculations and flow visualization experiments have shown that for 

small values of the Reynolds number. Re, transverse rolls occurred when the ratio of the 

Grashof number (Gr) to the Reynolds number, Gr/Re, was greater than a critical value that was 

T - T in 

a weak function of the dimensionless temperature difference, r T = — — ■= , where T sus is 

* in 

the temperature of the heated susceptor and is the temperature at the reactor inlet [2]. It was 
concluded that for Re < 10, Gr/Re should be kept well below 60 to prevent transverse rolls. 
Another study used two-dimensional simulations verified by experimental results to estimate 
the extent of transverse rolls (also referred to as recirculation) in a horizontal reactor [3]. This 
study determined that no transverse rolls should be expected for Gr/Re < 24. For 24 < Gr/Re 
< 75, the simulations showed areas of low velocity or nearly stagnant regions of gas that were 
replaced by transverse rolls in the same areas when Gr/Re > 75. Another computational study 

_3 

found that transverse rolls were always absent when Gr/Re < 100 for 1x10 < Re ^ 4 [1], 
The Grashof and Reynolds numbers were computed using fluid properties at the inlet 
conditions. These criteria yielded good agreement with reported experimental observations. 

Some researchers have su gg ested the use of Gr/Re 11 to determine the presence of transverse 
rolls, where n =1 for small Reynolds numbers and n=2 for larger Reynolds numbers. 
Transverse rolls were found to occur when Gr/Re 2 was greater than some critical value which 
depended on the dimensionless temperature difference, as defined above [2] . This criteria 
applied for Re £ 8 . In another study [1], transverse rolls were not present when Gr/Re < 25 
for 4 £ Re < 100 , with the Grashof and Reynolds numbers computed using fluid properties at 
the inlet conditions. However, a number of other studies [3-6] have found that Gr/Re is not 
reliable in describing the flow in a horizontal reactor. 

The Rayleigh number, Ra, may be interpreted as the ratio of the thermal flux by free 
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convection to that by diffusion. It characterizes the convective state of an enclosed gas 
subjected to a destabilizing temperature gradient, such as a heated bottom and cold upper wall 
[7]. This is often the situation in horizontal CVD reactors with water-cooled walls and a heated 
deposition surface. The Rayleigh number has been used to predict the onset of longitudinal 
rolls in the flow. Longitudinal rolls have been found to appear for Ra >1708 [2,8]. 

In this study, expressions for Gr/Re 11 , n = 1 and 2, were derived in terms of common reactor 
operation parameters. Values for these ratios were then computed for an example reactor for 
two carrier gases, argon and hydrogen, at pressures in the range 2-101 kPa, and predictions 
were made as to the pressure at which transverse rolls should and should not be present. A 
similar analysis was carried out using values of Ra, from an expression in terms of reactor 
operation parameters, to determine the pressures at which longitudinal rolls would and would 
not be present. These predictions were then compared to a numerical simulation of the flow in 
the sample reactor, with particular attention paid to the flow over the deposition surface. 


2. ANALYSIS OF SYSTEM 

While dimensionless numbers cannot provide detailed information regarding a process 
as complex as CVD, they can provide useful information on the relative magnitudes of various 
effects. Dimensionless numbers are usually expressed in terms of viscosity, thermal 
conductivity, heat capacity, and density, which are not the parameters usually controlled in a 
CVD reactor. Expressions for some dimensionless numbers useful in CVD analysis have been 
derived in terms of parameters that would normally be controlled in a CVD reactor [9]. These 
expressions were obtained by assuming ideal gas behavior (a good assumption at the relatively 
high temperatures and low pressures used for CVD), and estimation techniques for viscosity, 


thermal conductivity, and heat capacity based on kinetic theory. From this analysis, the 
Grashof number, which may be interpreted as the ratio of the momentum flux by free 
convection to that by diffusion, may be expressed as: 


Gr = 2.031xl0 3 g(MCT 4 £2 2 )(L 3 ) 


^ATP 2 ^ 


( 1 ) 


When the gravitational constant, g, is equal to 9.8 m/s, equation (1) becomes: 


Gr = 1.9903 x10 4 (Mct 4 Q 2 )(L 3 ) 


ATP 

T 


(2) 


where: M = molecular weight (kg/kgmoles) 

a = hard sphere diameter (angstroms) 

Q = collision integral (dimensionless) 

L = characteristic length (m) 

AT = characteristic temperature difference (K) 
P = pressure (Pa) 

T = temperature (K) 


The Reynolds number may be expressed as [9]: 


Re = 3.747x10' 


2 ~ 

ct_s2 


Q. W L ^ m 

0 5 UJ T 0i 


(3) 
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where: cr = collision diameter (angstroms) 

Q = collision integral (dimensionless) 
M = molecular weight (kg/kgmole) 

L = characteristic length (m) 

= cross-sectional area (m ) 
m = mass flow rate (kg/s) 

T = temperature (K) 


2 

By combining equations (2) and (3), an expressions for Gr/Re and Gr/Re may be 
obtained: 


^ = 5.312x10 2 (M L5 a 2 Q)(L 2 A ) 
Re 


^ atp 2> 

mT 35 ; 


(4) 


Gr 



= 1.418x10 7 (M 2 )(LA 2 s ) 


^ atp a 

VmV> 


(5) 


where: M = molecular weight (kg/kgmole) 

ct = collision diameter (angstroms) 

Q. = collision integral (dimensionless) 

L = characteristic length (m) 

Acs = cross-sectional area (m ) 

AT = characteristic temperature difference (K) 
P = pressure (Pa) 
m = mass flow rate (kg/s) 

T = temperature (K) 
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k B T 

The collision integral, H , is a function of the dimensionless temperature, , as shown 

in Figure 2 [10]: 

The Rayleigh number, which may be interpreted as the ratio of the thermal flux by free 
convection to that by diffusion, characterizes the convective state of an enclosed gas subjected 
to a destabilizing temperature gradient, such as a heated bottom and cold upper wall [2, 7, 8]. 
This is often the situation in horizontal CVD reactors with water-cooled walls and a heated 
deposition surface. Longitudinal rolls have been found to appear for Ra >1708 [2]. The 
Rayleigh number for a monatomic gas may be expressed as [9]: 


Ra = 1.354xl0 3 g(MCT 4 Q 2 )(L 3 ) 


A ATp 2\ 


^ t 4 ; 


(monatomic gas) (6) 


When g = 9.8 m/s, equation (6) becomes: 


Ra = 1.3267 x10 4 (Mct 4 Q 2 )(L 3 ) 


^ATP^ 


(monatomic gas) 


V i ) 


( 7 ) 


The Rayleigh number for a diatomic gas may be expressed as [9]: 


Ra = 1.402xl0 3 g(Mo 4 Q 2 )(L 3 ) 


^ATP^ 


V 1 ) 


(diatomic gas) (8) 


When g = 9.8 m/s, equation (8) becomes: 


Ra = 1.3739xl0 4 (Ma 4 Q 2 )(L 3 ) ^ 


2 \ 


(diatomic gas) (9) 
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The constanst in equations (1) through (9) depend upon the units used for the other 
variables, so care must be taken when using these equations that all variables are in the 
specified units. 

As an estimate of the temperature in the reactor, the following equation is suggested: 

_ T$ub * ^wall (IQ) 

2 

where: T sub = substrate surface temperature (K) 

T wa u = reactor wall temperature (K) 

The temperature difference can be expressed as: 

AT = T aJb -T w-1 (ID 

where: T sub = substrate surface temperature (K) 

T W all = reactor wall temperature (K) 

Since the main area of interest in a CVD reactor is the deposition surface, equations (10) 
and (11) should be reasonable estimates of conditions above the surface. 

The entrance effects in horizontal reactors may be estimated from the following 
expressions [7]: 


x h = 0.04hRe 

(hydrodynamic entrance length) 

(12) 

Xt = 0.4hPrRe 

(thermal entrance length) 

(13) 
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where: x h = hydrodynamic entrance length (m) 

x T = thermal entrance length (m) 
h = channel height (m) 

Re = Reynolds number (dimensionless) 

Pr = Prandtl number (dimensionless) 

In this study, the characteristic length, L, has been chosen to be the channel height. 

Equations (12) and (13) then become: 


x h = 0.04LRe 

(hydrodynamic entrance length) 

(14) 

x T = 0.4LPrRe 

(thermal entrance length) 

(15) 


The Prandtl number may be expressed as: 



where: |I = Newtonian viscosity (kg/m s) 

Cp = constant pressure heat capacity (J/kg K) 
k = thermal conductivity (W/m K) 


(16) 


Using estimation methods for viscosity, heat capacity, and thermal conductivity [9], the 
Prandtl number is found to be approximately 0.67 for a monatomic gas and 0.69 for a diatomic 
gas. Equations (15) then becomes: 
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(17) 


x T = 0.268LRe (thermal entrance length - monatomic gas) 

x T = 0.276LRe (thermal entrance length - diatomic gas) (18) 

Equation (14) for the hydrodynamic entrance length applies to an isothermal channel 
of height L. This entrance length will be somewhat longer for the non-isothermal case used in 
this study [7]. Equation (15), and, therefore, equations (17) and (18), apply to a horizontal 
parallel-plate reactor heated from below [7], and should, therefore, be applicable to this study. 


3. MODEL 


In this study, a 0.515 m long reactor with a 0.0508 m square cross-section was used. A 
diagram of the modeled reactor is shown in Figure 1. The model has been described in detail 
previously [11]. Some changes were made to this model for this paper. Since the flow patterns 
alone were of interest for this study, no chemical reactions were included. Radiation heat 
transfer was added to the model. The susceptor temperature no longer remains constant, but 
varies. A constant volumetric heat rate of 1.53x10® W/m 3 was specified at the susceptor 
surface, resulting in a maximum temperature of 1823 K. The inlet volumetric flow rate was 
increased from 1.2xl0' 5 to 2.2x10-5 m 3 /s to more accurately match experimental conditions 
that have been used in this reactor. 

Axial and transverse velocities of zero were specified at the surface of the susceptor, the 
quartz support, and the reactor walls. The inlet and reactor walls were maintained at a constant 
temperature of 280 K. The temperature distribution in the quartz support was calculated using 
the thermophysical properties of quartz. Pressures ranging from 2 kPa to 101 kPa were 
modeled. The inlet molar flow rate was held fixed as the pressure was varied, resulting in an 
increase in the inlet velocity as the pressure was lowered. 


4. RESULTS AND DISCUSSION 


The properties of argon and hydrogen used for this analysis are shown in Table 1, with 
system conditions are shown in Table 2. The characteristic length was chosen to be the distance 
between the susceptor surface and the top reactor wall. Equation (14) was used to calculate the 
hydrodynamic entrance length for both argon and hydrogen carrier gases. Equation (17) was 
used to calculate the thermal entrance length for an argon carrier gas, and equation (18) was 
used to calculate the thermal entrance length for a hydrogen carrier gas. Since there seems to 
be some disagreement in the literature as to which ratio is appropriate for describing flow in a 
horizontal reactor, the Gr/Re and Gr/Re 2 values for both argon and hydrogen carrier gases were 
calculated to determine which ratio best described the presence or absence of transverse rolls 
in the reactor modeled for this work. The values computed for Gr/Re, Gr/Re 2 , and Ra are given 
in Table 3 for argon and Table 4 for hydrogen. Equation (4) was used to compute Gr/Re for 
both gases, equation (5) was used to compute Gr/Re 2 for both gases, equation (7) was used to 
compute Ra for argon, and equation (9) was used to compute Ra for hydrogen. The computer 
program FLUENT was used for the simulations [12], Velocity vectors were used to visualize 
transverse rolls, and these vectors were scaled differently for each case for clarity. Streaklines 
(massless particle trajectories) were used to visualize the longitudinal rolls in the reactor 
predicted by the computational model. Particles were introduced into the reactor at a height 
of four centimeters from the bottom of the reactor, which is just above the plane of the 
deposition surface, and the trail of these particles as they followed the flow was recorded. 

The hydrodynamic entrance length was found to be about 0.03 meters for an argon 
carrier gas and about 0.004 meters for a hydrogen carrier gas. The velocity profiles will be 
changing over a larger portion of the deposition surface when argon is used as a carrier gas 
compared to hydrogen. In the sample reactor used in the simulations, the velocity profiles will 
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be changing over about the first third of the susceptor. The thermal entrance length was found 
to be about 0.19 meters for an argon carrier gas and about 0.025 meters for a hydrogen carrier 
gas When an argon carrier gas is used, the temperature profiles will be changing over the entire 
length of the deposition surface examined in this study. The temperature profiles will only be 
changing over about the first third of the deposition surface when hydrogen is used as the 
carrier gas. With such a large thermal entrance length for argon, large axial temperature 
variations will likely occur over the deposition surface, raising the question as to whether 
dimensionless numbers can accurately predict flow patterns for this carrier gas. As discussed 
later, predictions from these numbers proved valid for argon. 

Since the inlet molar flow rate was held constant, the Reynolds number remained the 
same for each gas even as the pressure was varied. For an argon carrier gas, the Reynolds 
number was 30. For a hydrogen carrier gas, the Reynolds number was 4. The values for Gr/Re 
were plotted for both carrier gases in Figure 3. From these plots, the approximate pressures 
corresponding to critical values for the dimensionless numbers can be determined. Using the 
criteria from [3], no transverse rolls should exist for Gr/Re < 24, while transverse rolls should 
exist for Gr/Re > 75. From the plots of Gr/Re versus pressure in Figure 3a this implies that no 
transverse rolls should exist for an argon carrier gas at pressures less than about 26 kPa, while 
transverse rolls should exist at pressures above about 49 kPa. From Figure 3b, no transverse 
rolls should exist for a hydrogen carrier gas at pressures less than about 67 kPa. A value of Gr/ 
Re greater than 75 is not reached for hydrogen over the pressure range examined here. The 
computational simulation results, as shown by velocity vectors in Figures 5-10, agree with 
these predictions. For both carrier gases, transverse rolls begin to be seen at the next highest 
pressure examined above the pressure determined to be the limit below which transverse rolls 
should not occur. 

The ratio Gr/Re 2 did not yield useful predictions in this study. The values of this ratio 
at each pressure were the same for argon and hydrogen carrier gases. A plot of Gr/Re versus 
pressure is shown in Figure 4. From [1], no transverse rolls should exist when Gr/Re < 25 
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II 


when 4 £ Re < 100 . Both argon and hydrogen carrier gases had Reynolds numbers within this 
range at all pressures examined, although hydrogen’s value was right on the lower limit. Both 
argon and hydrogen have values of Gr/Re^ below 25 at all pressures examined, yet transverse 
rolls existed in the simulations. 

As mentioned previously, longitudinal rolls should exist for Ra > 1708. As seen from 
Figure 11, the pressure corresponding to this Rayleigh number is about 54 kPa for an argon 
carrier gas; the pressure for hydrogen corresponding to this critical Rayleigh value was not 
reached in the pressure range examined. The computational model results, as shown in Figures 
12-17, agree with these predictions, except for hydrogen at 101 kPa, where longitudinal rolls 
are clearly present. Some waviness in the flow above the deposition surface is observed for the 
argon carrier gas at a pressure of 40 kPa and for a hydrogen carrier gas at a pressure of 80 kPa, 
but well-developed longitudinal rolls are still not present. 

There are errors associated with the approximations used to obtain equations (7) and (9) 
[9]; however, these errors are not large enough to result in the order of magnitude difference 
which would be necessary for the hydrogen carrier gas at 101 kPa to reach critical Rayleigh 
number. The critical Rayleigh number of 1708 has been obtained for a Boussinesq fluid 
between parallel plates of constant temperature, the bottom plate at a higher temperature than 
the top one [8, 13, 14], The Boussinesq approximation assumes a constant-property fluid, 
except for density in the gravitational term of the momentum equation, where it assumes a 
simplified equation of state to approximate the density difference [15]. In the simulated CVD 
reactor, the fluid properties were variable with temperature, and the susceptor surface (the 
“bottom wall”) although at a higher temperature than the top wall, was not at a constant 
temperature. The reactor also had a finite length and sidewalls. All of these differences could 
have contributed to the difference between the predictions regarding longitudinal rolls at 101 
kPa for hydrogen and the simulation results. In addition, dimensionless numbers must be 
calculated from conditions at a particular location, despite the fact that these conditions can 
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change greatly throughout the reactor. Therefore, these numbers can only provide us an 
approximate picture of the complex conditions inside a CVD reactor. 


5. CONCLUSIONS 

The ratio of the Grashof number to the Reynolds number, Gr/Re, has been found to be 
useful for predicting the pressures at which transverse rolls will be present in a horizontal CVD 
reactor. Predictions made from the values of this ratio agreed with those from a numerical 
simulation of a CVD reactor for argon and hydrogen carrier gases over the pressure range 2- 
101 kPa. The ratio Gr/Re2 was not found to be useful in predicting the pressures at which 
transverse rolls will be present in a horizontal CVD reactor. 

The Rayleigh number, Ra, was found to be useful for predicting the pressure at which 
longitudinal rolls will be present in a horizontal CVD reactor. Predictions made from the 
values of this number agreed with numerical simulations of a horizontal CVD reactor for the 
carrier gases argon and hydrogen over the pressure range 2-101 kPa. The only exception was 
a hydrogen carrier gas at 101 kPa. The calculated value of the Rayleigh number predicted that 
longitudinal rolls would not be present; however, the numerical simulation clearly showed that 
such rolls were present. It must be remembered that the critical value of 1708 for the Rayleigh 
number used in predicting the presence of longitudinal rolls was derived for a set of conditions 
that are usually not met in a CVD reactor. This could at least partially explain the discrepancy 
between the prediction from the Rayleigh number and the numerical simulations at 101 kPa. 
Also, dim ensionless numbers are calculated from conditions at a specific location, despite the 
fact that these conditions can vary greatly throughout the reactor. In such situations, 
dimensionless numbers cannot provide an exact description of a CVD reactor, and care must 
be used in interpreting the predictions obtained from them. 
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NOTATION 


Symbol 

Description 

Units 

Acs 

cross-sectional area 

_2 

m 

C P 

constant pressure heat capacity 

J/kg K 

Dim 

diffusivity of species i in a mixture 

m 2 /s 

8 

gravitational constant 

m/s 2 

h 

channel height 

m 

k 

thermal conductivity 

W/mK 

k B 

Boltzmann constant 

J/K 

Gr 

Grashof number 

dimensionless 

L 

characteristic length 

m 

m 

mass flow rate 

kg/s 

M 

molecular weight 

kg/kgmoles 

P 

pressure 

Pa 

Pr 

Prandtl number 

dimensionless 

r T 

dimensionless temperature difference 

dimensionless 

Ra 

Rayleigh number 

dimensionless 

Re 

Reynolds number 

dimensionless 

Sc 

Schmidt number 

dimensionless 

T 

temperature 

K 

T su b 

substrate surface temperature 

K 

T wa ii 

reactor wall temperature 

K 

AT 

characteristic temperature difference 

K 

x h 

hydrodynamic entrance length 

m 

*t 

thermal entrance length 

m 
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a 

e 

a 

v 

Q 


2 

thermal diffusivity m /s 

characteristic energy of interaction between molecules J 

Newtonian viscosity kg/m s 

hard sphere diameter angstroms 

kinematic viscosity m / s 

collision integral dimensionless 
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Table 1: Properties of argon and hydrogen at T = 1051.5 K [16] 


Parameter 

Argon Value 

Hydrogen Value 

M 

39.944 

2.016 

a 

3.418 

2.915 

e/k 

124 

38 

Q 

0.846 

0.710 


Units 

kg/kgmoles 

angstroms 

K 

dimensionless 
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Table 2: System Conditions Used for Analysis 


Parameter 

L 

Acs 
A T 
T 

HI argon 
m hydrogen 


Value 

Units 

0.0218 

m 

l.llxlO' 3 

m 2 

m 

1543 

K 

1051.5 

K 

8.693xl0' 5 

kg/s 

4.387xl0’ 6 

kg/s 


20 


ti 



Table 3: Values for Gr/Re, Gr/Re 2 , and Rayleigh numbers at various pressures: argon carrier 
gas (Re = 30) 


Pressure (kPa) Gr/Re 

2 lxlO' 1 

20 lxlO +1 

40 5xl0 +1 

60 lxlO +2 

80 2xl0 +2 

101 3xl0 +2 


Gr/Re 2 

Ra 

4xl0' 3 

3x10+° 

£ 

o 

1 

3xl0 +2 

2xlO +0 

lxl0 +3 

4xlO +0 

2xl0 +3 

7x10+° 

4xl0 +3 

lxl0 +1 

7xl0 +3 
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Table 4: Values for Gr/Re, Gr/Re 2 , and Rayleigh numbers at various pressures: hydrogen 
carrier gas (Re = 4) 


Pressure (kPa) 

Gr/Re 

Gr/Re 2 

Ra 

2 

2xl0' 2 

4xl0' 3 

5xl0' 2 

20 

2x1 0 +0 

4X10" 1 

5x10+° 

40 

7x10+° 

2x10+° 

2xl0 +1 

60 

2xl0 +1 

4xlO +0 

5xl0 +1 

80 

3xl0 +1 

7x10+° 

8xl0 +1 

101 

5xl0 +1 

lxl0 +1 

lxl0 +2 
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Figure 1: CVD reactor used in model. All dimensions in meters. 


Figure 2: Collision integral for viscosity and thermal conductivity calculations as a function of 

k B T k B T 

dimensionless temperature for the range: a) 0 < < 10; b) 10 < — — < 400. 


Figure 3: Gr/Re versus pressure for: a) argon; b) hydrogen 

Figure 4: Gr/Re^ versus pressure for argon and hydrogen. The values of this ratio were the 
same for both gases. 

Figure 5: Velocity vectors showing transverse rolls at 2 kPa for : a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 6: Velocity vectors showing transverse rolls at 20 kPa for: a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 7: Velocity vectors showing transverse rolls at 40 kPa for: a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 8: Velocity vectors showing transverse rolls at 60 kPa for: a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 9: Velocity vectors showing transverse rolls at 80 kPa for: a) argon; b) hydrogen. Flow 
enters from the left. 


Figure 10: Velocity vectors showing transverse rolls at 101 kPa for : a) argon; b) hydrogen. 
Flow enters from the left. 
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Figure 1 1: Ra versus pressure for: a) argon; b) hydrogen 


Figure 12: Streaklines showing longitudinal rolls at 2 kPa for : a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 13: Streaklines showing longitudinal rolls at 20 kPa for: a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 14: Streaklines showing longitudinal rolls at 40 kPa for: a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 15: Streaklines showing longitudinal rolls at 60 kPa for: a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 16: Streaklines showing longitudinal rolls at 80 kPa for: a) argon; b) hydrogen. Flow 
enters from the left. 

Figure 17: Streaklines showing longitudinal rolls at 101 kPa for a) argon; b) hydrogen. Flow 


enters from the left. 
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